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Abstract 

We discuss a possibility that the so-called reactor antineutrino anomaly can 
be, at least in part, explained by applying a quantum held-theoretical ap¬ 
proach to neutrino oscillations, which in particular predicts a small deviation 
from the classical inverse-square law at short but macroscopic distances be¬ 
tween the neutrino source and detector. An extensive statistical analysis of 
the reactor data is performed to examine this speculation. 


Introduction 


Nuclear reactors are intense sources of electron antineutrinos whose spec¬ 
trum is composed of thousands of spectral components formed mainly by 
the P decay of the hssion products of the four parent isotoj^s: 

^^®Pu, and ^^^Pu. The very sophisticated recent calculations [l|-l3| yield a net 
3-3.5% upward shift in the predicted energy-averaged Ve flux with respect to 
the previously expected flux used in the earlier short baseline (^BL) reactor 
experiments (ILL-Grenoble jl, , Gosgen , Krasnoyarsk 0, [sj, Rovno ji- 
Tl| . Bugey [m, 1^, Savannah River Plant |14|1. The Uf, flux normalization 
uncertainty in the new calculations is claimed to be only ±2.7%. This im¬ 
plies 0 that the measured event rates in the SBL experiments are about 6% 
too low, giving rise to the so-called “reactor antineutrino anomaly” (RAA). 

Figured] illustrates this issue. The curve shows the ratio of the W-induced 
event rate calculated with and without regard for the oscillations. Here 
and thereafter we use the global best £t values for the neutrino mass-squared 
splittings and mixing angles from Ref. [lH for the normal mass hierarchy; 
we also assume no CP violation in mixing. The cross section for the in¬ 
verse d decay (IBD) is calculated by using the recent analytical results of 
Ref. 1^, which take into account the radiative corrections of order a/n and 
contributions of weak magnetism and neutron recoil to next-to-leading or¬ 
der in the expansion in inverse powers of the nucleon mass. In Fig. [T] and 


^Talk at the International Workshop on Prospects of Particle Physics: 
rino Physics and Astrophysics”, Valday, Russia, February 1-8, 2015. 
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Figure 1: The reactor antineutrino anomaly. The points represent the ratios of 
measured event rates to that expected with no oscillations. The vertical error 
bars do not include the common normalization uncertainty. The original data 
from Ref s. |^ 14l. I17 h 21I] are recalculated as explained in the text. The data from 
Refs. 22, 2^ are relative. The ratios of the KamLAND 21] and Daya Bay 2015 24] 
results to the “Huber-Mueller” model predictions are plotted at the flux-weighted 
distances of 175 km and 573 m, respectively. The curve represents the theoretical 
expectation, in which the Sv oscillations parameters are fixed to the global best fit 
values taken from the Review of Particle Physics 15|]; the normal mass ordering 
and no CP violation are assumed. 


below, all the curves correspond to a reactor with pure fuell The 
original data from Refs. js-Jd, 17-21] are corrected according to Ref. [l| 
and then renormalized to the new world average value of the neutron mean 
life [HI. It is seen that most of the SBL data points (the measurements at 
L < 100 m, where L is the distance between the reactor core and detec¬ 
tor) are below the expectation. A clear trend is visible at L < 20 m that the 
closer detector is located to reactor, the smaller the measured rate (the larger 
the discrepancy between the data and theory). Note that the data points 
“Krasnoyarsk-IV” [^ and “Rovno 92” [ll| are ignored in the numerous RAA 
analyses, but the latter point is signihcant for revealing the mentioned trend. 


^Nevertheless, in the following analysis we explicitly take into account the particular 
fuel composition in each experiment, since it does have a small impact on the results. 
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As is seen, the theory is in rather poor agreement with the latter result. 
Hence, both the earlier SBL and new Daya Bay measurements give a hint to 
either “new physics”, or merely a lower Vf. flux than predicted in Refs. UM- 


2 Extra neutrinos or wrong normalization? 


Most if not all efforts to resolve the anomaly are based on the hypothesis 
of existence of one or more light (eV mass scale) sterile neutrinos, that is 
fundamental neutral fermions with no standard model interactions except 
those induced by mixing with the standard (active) neutrinos. The active-to- 
sterile neutrino mixing would lead to a distance-dependent spectral distortion 
and overall reduction of the reactor Vf> flux. 

In Fig. |2]we show, as an example, the results of calculations performed 
in the framework of the simplest “3-1-1” phenomenological model with one 
sterile (anti)neutrino, 1 ^ 4 , by using the three pairs of the z /4 — ui mixing 
parameters, (Am 4 i, sin^ 26 * 41 ), listed in the legend of the figure. These values 
were derived in Ref. from detailed statistical analyses of all the neutrino 
oscillation data available to date. The “SBL rates only” fit, includes the SBL 
reactor data except the points “Krasnoyarsk-IV” and “Rovno92” (see Fig. [T]). 
The “SBL -I- BugeyS spectrum” fit, includes the same data set and spectral 
data from BugeyS 13|. The “Global Uf, di sapp earance” fit involves the data 
from the reactor experiments jl,[6-10, 12-14, 17-19, |^, as well as 

solar neutrinos (261 data points from Homestake, SAGE, GALLEX/GNO, 
Super-Kamiokande, and SNO experiments), radioactive source experiments 
at SAGE and GALLEX, and the LSND/KARMEN Ue disappearance data 
from — scattering (see Ref. for the full list of references and further 
details). It is necessary to note that these fits operate with somewhat lower 
(to within roughly 1 %) values for the reactor 17^ induced rates and with a 
little bit different covariance matrix, as compared to those used in the present 
analysis (the details of our calculations will be published elsewhere). 

The solid curve in Fig. [2] represents the same 3z/ oscillation prediction 
as in Fig. [H but shifted down by the normalization factor Nq derived from 
a fit to all the data except these from RENO and Daya Bay. In this fit 
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Figure 2: Comparison of the data with the “3+1” model (see text) and with the 
standard 3i^ oscillation prediction after renormalization of the Vf. flux (solid curve). 
The renormalization factor is obtained from a fit to all the data but RENO and 
Daya Bay. The filled band shows the ±1(7 uncertainty of the fit. The data points 
are the same as in Fig. [T] but the eight relative data points are shifter by the factor 
iVo = 0.950. 


we take into account the correlation between the data, including the over¬ 
all normalization uncertainty, which is taken to be 2.7% j^. The obtained 
factor TVo = 0.950 ^q;o26 (x^/NDF 1) does not contradict to the adopted 
flux uncerta inty , but is somewhat different from the results of earlier calcula¬ 
tions j^, 16, 271, which used different data subsets and input parameters. All 
four curves in Fig. p are in agreement, within the errors, with the new Day 
Bay measurement [2^, but are in some conflict with the ILL data point j^. 

Recently, it was argued that the true uncertainty in the Fg flux pre¬ 
dictions may be as large as 5% and the spectral shape uncertainties may be 
much larger due to poorly known structure of the forbidden decays. This 
hnding has been in essence conhrmed by the new precision measurements 
of the positron energy spectra from IBD which show apparent ~ 10% 
excess in 4 — 6 MeV region compared to the expectation based on the models 
of Refs. [3,3 (see Ref. for further discussion and references). From what 
has been said it appears that the efforts to explain the anomaly by the sterile 
neutrino hypothesis may be somewhat premature. Moreover, it is seen from 
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Fig. |2] that the proper renormalization of the flnx is hardly distingnishable 
from the “global v^. disappearance” £t and (maybe somewhat accidentally) 
almost fnlly coincides with the “SBL + BngeyS spectrnm” ht. We empha¬ 
size however that the steady decrease of the event rate at very small L, ij 
real, cannot be explained by neither the wrong W flnx normalization alone, 
nor the “3-1-1” scenario. Thns, it is appropriate to consider an alternative 


explanation. Snch an alternative has been proposed in Ref. jSOj. It is based 


on a qnantnm field-theoretical (QFT) approach to the nentrino oscillation 
phenomenon, which predicts a small deviation of the (anti)nentrino event 
rate, as a fnnction of the distance L between the sonrce and detector, from 
the classical inverse-sqnare law (ISL) behavior. Below we consider this issne 
in some detail. 


3 A sketch of the QFT approach 

The “nentrino-oscillation” phenomenon in the F-matrix QFT approach is 
nothing else than a resnlt of interference of the macroscopic Feynman di¬ 
agrams (like shown in Fig. [3]) which describe the lepton nnmber violating 

Fignre 3: A generic macroscopic Feyn¬ 
man diagram. Here Is,d and Fs^d de¬ 
note the sets of the initial (/) and final 
(F) WP states in the “source” (W) 
and “detector” {Xd) vertices; F'^ ^ = 

; where are the charged 
lepton WP states (a,/3 = e,fi,T); qs^d 
are the 4-momentum transfers in the 
vertices, as defined by Eq. ([5]). The 
vertices are in general macroscopically 
separated in space and time. The par¬ 
ticular “decryption” of the neutrino 
production/absorption mechanism as¬ 
sumes the standard model charged cnr- 
rent interaction of quarks and leptons. 

processes with the neutrino mass eigenhelds z/j (i = 1,2,3) treated as inter¬ 
nal lines (propagators). The external lines of the macrodiagrams are assumed 
corresponding to asymptotically free quasi-stable wave packets (WP) rather 
than the conventional to QFT one-particle Fock’s states |k, s) with definite 
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3-momenta k and spin projections s. According to Refs. 3l|, l32|, the free ex¬ 


ternal WP states are constructed as covariant space-time point x dependent 
linear superpositions of the one-particle states, 


|p,s,a:) = j 


(27r)32P;k 


k,s), 


( 1 ) 


satisfying the correspondence prineiple which demands that |p, s,a:) turns 
into |k, s) in the plane-wave limit (PWL) that is equivalent to the following 
condition for the relativistic invariant form factor function 0: 


0(k,p)^(27r) 


^2E^6{k - p) 


The detail properties of the WP states ([I]) are discussed in Refs. [33|, [SJ . 

Within the outlined approach, and after applying several more or less 
“technical” simplihcations, it is proved that the neutrino induced event 
rate in an ideal detector can be written (somewhat symbolically) in the form 


^ J dx J dy j d^y j da^vVafs |y - x|). (2) 

Here r is the detector exposure time, E^, is the neutrino energy, Va/s is 
the QFT generalization of the standard quantum-mechanical neutrino flavor 
transition probability, the differential form da^D represents the differential 
cross section of the neutrino scattering from the whole detector device; d^v 
is the differential neutrino flux incident on the detector from a stationary 
source device (e.g., a hssion reactor core). The integrations in ([2]) are over 
the source and detector hducial volumes Vs and Vt>. The theory explicitly 
predicts that the neutrino flux decreases with increasing the distance |y — x| 
in compliance with the usual inverse-square law (ISL): 


ddu oc |y - X 


(3) 


This quite expected result has been derived by using the so-called Grimus- 
Stockinger (GS) theorem 3^, which dehnes the asymptotic behavior of the 


amplitude at L —?■ oo and this is the crucial point in the context of the 
problem under consideration. 

As it follows from the formalism, the L dependence of the amplitude 
described by the macrodiagram shown in Fig. [3] is dehned by the neutrino 
propagator modified by the external wave packets. 


id + m, 


d*q S,(q - q,)Sd(q + gd)e 


(2ir) 


— rnf + ie 


(4) 
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where x = (|/o — Xq, y — x) = (T, L), and are the 4-momentum transfers, 
(ls = ^Pa-^Pb, qd=^Pa-^Pb, (5) 

q-G/s b^Fs CL^Id b^Fd 

are the most probable (on-shell) 4-momenta of the external packets x G 
Is ® Id ® Fg ® Fdi and mj is the mass of the neutrino held z/j. The functions 
(5s(g — Qs) and 5d{q + qd) are the “smeared” 5 functions (see Ref. for 
their explicit form) dehned by the 4-momenta masses {rril^ = p^'^ 


k) 5 
2 ^ 


and momentum spreads (T^ of the external in and out packets {a\ <g; m 
In the plain-wave limit ((T^ —?■ 0, Vx) these functions turn into the ordinary 

Dirac (5 functions, Isiq - qs) (5(g - qs), Sd{q + qd) S{q + qd), thus 
leading to the exact energy-momentum conservation in the vertices of the 
macrodiagram, and the function (jH) becomes, up to a multiplier, the standard 
fermion propagator. If however the momentum spreads are hnite, the 
space-time behavior of the function (jl]) is nontrivial. In particular, its spatial 
dependence at sufficiently large distances L is given by the above-mentioned 
GS theorem 3^, according to whichU 


J(L,k) = 


as L = |L| 


(iq <h(q)e' 


zqL 




$(-k1) 


(2vr) 


q2 _ ^2 _ 


4:71L 


l + O 


Vl 




—)■ cxD. This offers the QFT explanation of the ISL behavior ([3]) 
but does not, however, provide the spatial scale above which the distance L 
may be considered as “sufficiently large”. 

In Ref. [^, an extended version of the GS theorem has been proved, 
which parametrically dehnes such a scale by using the asymptotic expansion 
of the integral J(L, k) in terms of inverse powers of L at large L. To be more 
precise, the theorem in its simplest form states that for any function <h(q) in 
the Schwartz space S'(R^) 


J{L,k) = 


Ak,L 


AnL 


®(q) + 


(-i)^T>„<h(q) 


n>l 


L —y cxo, 


( 6 ) 


q=— k1 


where Dn are explicitly dehned differential operators on the momentum 
space; the lowest order operators, sufficient for our present purpose, are 

= I [V= - (1VJ=] - (IVJ, 

02=^ [V= - (IV" - k( 1V,) [V^ - - i [V= - 3(1V,)=] . 


^It is assumed that the complex-valued function ^(q) itself and its first and second 
derivatives decrease at least like l/|qp as |q| —>■ oo and > 0; see Ref. for details. 
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Additional important featnres can be fonnd in Refs. jSO, l3g. An analysis 
of Eq. ([n]) shows that the 1/L behavior of the amplitnde (and thns the ISL 
behavior of the event rate) is violated at the distances L < 2,0, where 


^eS 


(n^) (ft) ■ - 

and the fnnction aes = (JeffiK-, represents an effective momen- 

tnm spread dependent on the nentrino momentnm k as well as on the mean 
velocities, masses, and momentnm spreads of the external (in and ont) wave 
packets x. The explicit form of this fnnction can be fonnd after specihcation 
of a particnlar model for the external WP states. A simple example is dis- 
cnssed in Ref. l30l w ithin the so-called contracted relativistic Gaussian packet 
(CRGP) model |3,i. . It is in particular shown that aes is dehned through 

the transverse (with respect to the neutrino propagation direction 1 ) compo¬ 
nents of the inverse overlap tensors which determine the effective space-time 
overlap volumes of the WP states in the vertices of the macrodiagram. It 
is signihcant that these components are nearly independent of the neutrino 
masses (assuming these to be small with respect to the neutrino energy and 
thus K ~ Ey). Within the GRGP model, it can be also shown that the mag¬ 
nitude of (Jeff is strongly affected by the hierarchy of the external momentum 
spreads but in the simplest case when these spreads are similar in order 
of magnitude, aes is of the same order, too. So, as is seen from Eq. ([7]), the 
spatial scale ([7]) can be macroscopically large at sufficiently small external 
momentum spreads, thus leading to a measurable ISL violation (ISLV). 

It is shown in Ref. jssj that Eq. (|6]) modihes the formula for the event 
rate ([2]) in such a way that the relation ([3]) for the flux is replaced by 


^eff 


-2 


ddu OC 




1 2n 


( 8 ) 


where the coefficient functions are explicitly dehned from Eq. ([6]). By 
making expedient assumptions, it can be proved (and this is a crucial point) 
that Cl < 0. Hence, using Eq. ([H]) in leading order (thereby assuming that 
the ISLV correction is small), yields the following simple replacement for the 
event rate: 


dN f dN 

dr y j dr 


(9) 


_2 _ 

(provided that IIq L^). Here 2o ~ (£ 0 ) is a neutrino energy dependent 
parameter of dimension of length. Needless to say, at present this parameter 
cannot be obtained from hrst-principle calculations, but it can be measured. 







4 Data analysis 


To check the assumption that the ISLV effect could actually be, in part, 
responsible for RAA, we performed a statistical analysis of the reactor data 
discussed in Sect. [T] Since in this paper we use only the spectrum-averaged 
event rates, the L dependent factor in Eq. ([2]) can be replaced by (1 —Lq/L^), 
where Lq ~ (£ 0 ) is an energy independent parameter, which is a subject of 
the present study. Taking into account the large uncertainty in the flux 
normalization, we shall use the following theoretical model to fit the data: 


T{L;No,Lo) 


„ q _ Jo” St EME.)S,(E,) 

"V ib Jr‘iE.T.kMEu)Sk(E.) 


Here is the required normalization parameter, is the reactor fissile iso¬ 
tope fraction, Sk{E^) is the Rg energy spectrum (taken from Ref. [H), <j{Ey) 
is the IBD cross section [l^ and E^) is the Rg survival probability in 

the standard 3z/ mixing schemeEl In order to find the best-fit parameters Nq 
and Lq we minimize the standard with the full covariance matrix for the 
correlated data. The results of the fits of several types are shown in Figs. HHHl 



We thereby neglect the decoherence effects predicted in the QFT approach , reason¬ 
ably assuming that the baselines under consideration are too short for their manifestation. 
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Figure 5: Two-parameter fits to SBL data subset (a) and to all data (b). The 
data of RENO and Daya Bay are not included into the analysis. The eight relative 
data points in the panels (a) and (b) are shifter by the factor Nq = 0.975 and 
0.966, respectively. The inserts clarihes the agreement between the best-fit curves 
and Daya Bay data points. 
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Figure 6: Panel (a): The 68% C.L. error contours for the pair of the parameters Nq 
and Lq extracted from the fits to all data and SBL data subset (L < 300 m). The 
one-parameter (with Nq fixed to 1) fit to SBL data is also shown for comparison. 
Panel (b): The same as in (a), but with the ILL data point excluded from the 
statistical analysis. In both panels, the circles correspond to the best-fit values of 
the fitted parameters with and 

Figure m represents the results of the simplest one-parameter ht, in which 
Nq is set to 1 and only the SBL data (L < 300 m) are used in the analysis. 
Here and below, the hlled band represents the Tier uncertainty. It is seen 
that despite relatively small value of y^/NDF the best-£t curve is in rather 
poor agreement with the data; in particular, it is in conflict with the recent 
Daya Bay measurement. This indicates that the ISLV suppression alone is 
not sufficient and the ffux renormalization is actually required. 

Figure |5] illustrates the results of the two-parameter hts performed with 
two different data subsets, namely with the SBL data only, and with all the 
data except these from RENO and Daya Bay. As is seen from the figure, 
these two hts are in reasonable agreement to each other and both describe 
the reactor data rather well. Although the value of y^/NDF for the SBL 
data is nominally a bit better than that for the ht to all data (0.63 and 0.74, 
respectively), the latter ht is (as is clearly seen in the inserts of Fig. |5]) in 
better agreement with the Daya Bay point (let us remind that it does not 
participate in the analysis). By comparing these results with the hts shown in 
Fig. [21 we may conclude that the ISLV effect in combination with the proper 
renormalization of the Fg hux provides a better resolution of the anomaly. 
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To gain a deeper understanding of our results, we compare in Fig. E] 
the 68% C.L. error contours for the pair of the htted parameters (iVo,Lo), 
obtained from the hts to different data subsets. Panel (a) in Fig. [6] shows 
the contours for the hts to all data and SBL data subset. Panel (b) shows 
the same but with the ILL data point excluded from the analysis. It is seen 
with no need of additional explanation that the Grenoble experiment is the 
“cornerstone” for verihcation of the ISLV effect. It is however important that 
even without the ILL point we obtain essentially the same best-ht values of 
the parameters. 


5 Conclusions 


The QFT approach predicts a deviation from the classical inverse-square law 
at short baselines. While the numerical value of the spatial scale at which 
the deviation becomes essential cannot be predicted from the present-day 
theory, it can be extracted, under reasonable assumptions, from the data 
of the past and current reactor antineutrino experiments. Our statistical 
analysis demonstrates that the averaged over the reactor antineutrino spec¬ 
trum value of the scale (Lq) is about 3 m that roughly corresponds to the 
spectrum-averaged effective momentum spread ((Jeff) of about 0.5 — 0.8 eV 
(thereby hinting that the wave packets of the particles and nuclei involved 
into the reactor W production and detection may have “mesoscopic” effective 
dimensions). This is in agreement with the conservative estimate presented 
in Ref. js^. Besides, the best-ht value of Lq is very stable with respect to 
choice of the data subset and z/g spectrum model. To check the latter, we 
performed the same one- and two-parameter hts as described above, but with 
the input W energy spectra from Refs. ^ and 38|] ^rjwed by very diherent 

and the cumulative 


methods, as well as with combinations of the models 
Ve spectrum from gggjon recently measured with the scientihc neutron 
source FRMII in Garching . We conclude from these exercises that value 
of Lq is almost insensitive (within the errors) to the spectrum variations. 
Needless to say, it is not the case for the normalization parameter Nq. 

Although the available reactor data cannot dehnitely conhrm or exclude 
the light sterile neutrino hypothesis, and do not provide unambiguous sup¬ 
port for the ISLV ehect, they are in much better agreement with the latter. 
The next-generation experiments with very short baselines (L ;< 20 m), small 
neutrino or antineutrino sources, and high-precision, desirably movable de¬ 
tectors are required in order to conhrm or disconhrm our explanation. 
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